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ABSTRACT 

The discovery of transiting "super-Earths" with inflated radii and known masses such 
as Kepler-llb-f, GJ 1214b and 55 Cnc e, indicates that these exoplanets did not 
lose their nebula-captured hydrogen-rich, degassed or impact-delivered protoatmo- 
spheres by atmospheric escape processes. Because hydrodynamic blow-off of atmo- 
spheric hydrogen atoms is the most efhcient atmospheric escape process we apply a 
time-dependent numerical algorithm which is able to solve the system of 1-D fluid 
equations for mass, momentum, and energy conservation to investigate the criteria 
under which "super-Earths" with hydrogen-dominated upper atmospheres can experi- 
ence hydrodynamic expansion by heating of the stellar XUV (soft X-rays and extreme 
ultraviolet) radiation and thermal escape via blow-off. Depending on orbit location, 
XUV flux, heating efficiency and the planet's mean density our results indicate that 
the upper atmospheres of all "super-Earths" can expand to large distances, so that 
besides of Kepler-llc all of them experience atmospheric mass-loss due to Roche lobe 
overflow. The atmospheric mass-loss of the studied "super-Earths" is one to two or- 
ders of magnitude lower compared to that of "hot Jupiters" such as HD 209458b, so 
that one can expect that these exoplanets cannot lose their hydrogen-envelopes during 
their remaining lifetimes. 

Key words: planets and satellites: atmospheres - planets and satellites: physical 
evolution - ultraviolet: planetary systems - stars: ultraviolet - hydrodynamics 



1 INTRODUCTION 

The detection of hydrogen- and volatile-rich exoplanets at 
orbital distances <1 AU opens questions regarding their up- 
per atmosphere structures and the stability against escape 
of atmospheric gases. Since Js40% of all discovered exoplan- 
ets are orbiting their host stars at distances closer than the 
orbit of Mercury, the atmospheres of these bodies evolve in 
much more extreme environments than what is known from 
the planets in our Solar System. More intense stellar X-ray, 
soft X-ray, EU\|3 radiation and particle fluxes at such close 
orbital distances will alter the upper atmosphere structure 
of these objects to a great extent. 

Lammer et al. (2003) were the first to provide a model 
of exoplanetary upper atmospheres, which is based on an 



^ The radiation between 5-920 A contains soft X-rays and EUV 
and is considered as XUV radiation 



approximate solution of the heat balance equation in plan- 
etary thermospheres (Bauer 1971; Gross 1972) and found 
that hydrogen-rich upper atmospheres of Jupiter-type gas 
giants in close orbital distances will be heated to several 
thousands of Kelvin, so that hydrostatic conditions cannot 
be valid anymore, because the planet's upper atmosphere 
will expand dynamically upwards. Under such conditions 
the exobase can move to locations which are above a possible 
magnetopause or even at the Roche lobe distance (e.g., Lam- 
mer et al. 2003; Lecavelier des Etangs et al. 2004; Erkaev et 
al. 2007; Lammer et al. 2009a), so that the upward-flowing 
neutral gas can interact with the dense stellar plasma flow 
(Holmstrom et al. 2008; Ekenback et al. 2010; Lammer et 
al. 2011). Hydrodynamic models by Yelle (2004), Tian et al. 
(2005a), Garcia Munoz (2007), Penz et al. (2008), and Kosk- 
inen et al. (2010) agreed also with the hypothesis of Lammer 
et al. (2003; 2009a) that close-in hydrogen-rich exoplanets 
experience dynamic atmospheric expansion and outfiow up 
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to their Roche lobes with mass-loss rates in the order of 
~ (1-5) X 10^° gs"\ 

The first observational evidence of an XUV-heated ex- 
tended non-hydrostatic upper atmosphere was obtained for 
the hot gas giant HD 209458b, which orbits a Sun-like star 
at 0.047 AU (Vidal-Madjar et al. 2003; Ben-Jaffel 2007; 
Ben-Jaffel & Sona Hosseini 2010). Excess absorption in the 
Lyman-Q line of ~9-15 % during the transit was interpreted 
as being caused by a population of atomic hydrogen located 
beyond the Roche lobe, indicating a hot and extended up- 
per atmosphere. Subsequent detection of O I and C II by 
Vidal-Madjar et al. (2004), as well as Si III by Linsky et al. 
(2010), supports the hypothesis that these heavy elements 
are likely to be carried there from deeper atmospheric levels 
with the dynamical hydrogen flow. Metals have been found 
around the Roche lobe of the highly irradiated "hot Jupiter" 
WASP-12b, an inflated gas giant at ~0.023 AU, by Fossati 
et al. (2010). As heavier elements can only be dragged to 
the upper atmosphere regions along with upward flowing 
hydrogen atoms, these observations of heavy elements at 
large distances from the planet most likely confirm scenar- 
ios which were suggested by Sekiya, Nakazawa, & Hayashi 
(1980), Sekiya, Hayashi, & Nakazawa (1981), Hunten, Pepin, 
& Walker (1987) and Pepin (1991; 2000) that hydrodynamic 
escape from early Solar System planets like Venus could have 
generated Ne and Ar isotope ratios close to the observed val- 
ues in its present-time atmosphere and noble gas ratios sim- 
ilar to those derived for Earth's initial atmosphere. In view 
of atmospheric isotope studies it is generally accepted that 
hydrodynamic outflow driven by the XUV-active young Sun 
could have been responsible for the observed heavy isotope 
enrichment in the atmospheres of Venus, Earth and Mars 
(e.g., Lammer et al. 2008; and references therein). 

Therefore, depending on the activity of a planet's host 
star, its orbit location, mass and size, one can expect that 
some exoplanets experience high atmospheric mass-loss dur- 
ing their lifetimes (Lammer et al. 2003; Vidal-Madjar et al. 
2003; 2004; Lecavelier des Etangs et al. 2004; 2010; Lecave- 
lier des Etangs 2007; Yelle 2004; Baraffe et al. 2004; Tian 
et al. 2005a; Erkaev et al. 2007; Hubbard et al. 2007; Kho- 
dachenko et al. 2007; Garcia Muhoz 2007; Koskinen, Ayl- 
ward, & Miller 2007; Penz et al. 2008; Penz, Micela, & 
Lammer 2008; Penz & Micela 2008; Davis & Wheatley 2009; 
Lammer et al. 2009a; Murray-Clay, Chiang, & Murray 2009; 
Linsky et al. 2010; Fossati et al. 2010; Leitzinger et al. 2011; 
Owen & Jackson 2012). 

In this context, the discoveries of very close-in low- 
mass exoplanets like CoRoT-7b (e.g. Leger et al. 2009), 
and Kepler-lOb (e.g., Batalha et al. 2011) at < 0.02 AU 
raise the question whether these objects have always been 
so-called "super-Earths", or if are they just remnants of 
once more massive gaseous planets which lost their whole 
hydrogen envelopes (e.g., Valencia et al. 2010; Jackson et 
al. 2010; Leitzinger et al. 2011). However, besides rocky 
"super-Earths" which orbit at very close orbital distances 
such as CoRoT-7b, Kepler-lOb, Kepler-18b (Cochran et al. 

2011) and planets in the Kepler-20 system (Fressin et al. 

2012) , the Kepler space observatory discovered several low- 
density "super-Earths" closely packed in the Kepler- 11 sys- 
tem (Lissauer et al. 2011), whose mean densities indicate 
that they have rocky cores which are surrounded by signif- 
icant amounts of hydrogen/helium envelopes (e.g., Lissauer 



et al. 2011; Lammer et al. 2012; Lammer 2012; Ikoma & 
Hori 2012). Moreover, the mean densities of other transiting 
"super-Earths" , such as GJ 1214b (Charbonneau et al. 2009) 
or 55 Cnc e (Demory et al. 2012; Endl et al. 2012; Gillon et 
al. 2012), also indicate substantial envelopes of light gases 
such as hydrogen and He or H2O. For GJ 1214b, the current 
status of observational evidence supports mainly two pos- 
sible atmospheres, namely a H/He-dominated atmosphere 
with clouds and low methane content, or a H20-dominated 
atmosphere (Bean et al. 2010, 2011; Croll et al. 2011; Desert 
et al. 2011; de Mooij et al. 2012). For the large "super-Earth" 
55 Cnc e three different composition-based hypotheses can 
be found in the literature, within the first two possibilities 
are a rocky planet with a H/He dominated atmosphere, or a 
supercritical water planet (Demory et al. 2012; Gillon et al. 
2012). The third hypothesis is based on a recent study by 
Madhusudhan et al. (2012) where these authors suggest that 
the mass and radius of 55 Cnc e can also be explained by 
a carbon-rich solid interior made of Fe, C, SiC, and/or sili- 
cates and without a volatile envelope. Although, there may 
be this possibility of an alternative more exotic explanation 
for 55 Cnc e, the discoveries of these low-density hydrogen 
and/or H20-rich planets can be seen as an indication that 
many "super-Earths" even located close to their host stars 
will not lose their initial volatile-rich protoatmospheres dur- 
ing their lifetimes (Ikoma & Hori 2012; Lammer 2012). If 
so, the consequences of these findings are very relevant for 
the probability and evolution of Earth-type class I habitats 
(Lammer et al. 2009b; Lammer 2012) and habitability in 
general (Pierrehumbert & Gaidos 2011; Wordsworth 2012). 
Understanding the efficiency of atmospheric loss processes 
at hydrogen-rich "super-Earths" is therefore crucial. 

The aim of this study is to investigate the environmental 
conditions under which "super-Earths" with hydrogen-rich 
upper atmospheres undergo hydrodynamic blow-off. For this 
study we focus only on "super-Earths" for which the sizes 
and masses are more or less well determined. For that rea- 
son we do not include transiting "super-Earths" such as HD 
97658b, Kepler-9d, or small terrestrial planets such as those 
discovered in the Kepler-42 system with unknown or un- 
certain masses, or "super-Earths" with known masses but 
unknown radii. 

In Sect. 2 we describe the adopted stellar and planetary 
input parameters and the hydrodynamic upper atmosphere 
model which is applied for the investigation if Kepler- lib, 
Kepler-llc, Kepler-Ud, Kepler-lle, Kepler-Uf, GJ 1214b 
and experience atmospheric blow-off or less efficient Jeans 
escape. In the controversial case regarding the composition 
of 55 Cnc e, we assume in this study that the planet is 
surrounded by a hydrogen envelope, which would be the 
case if rocky core of the planet is surrounded by a H/He, or 
H2O dominated atmosphere as suggested by Demory et al. 
(2012) or Gillon et al. (2012). 

In Sect. 3 we apply our hydrodynamic model to these 
low-density "super-Earths" and study the response of their 
upper atmospheres to the stellar XUV flux and calculate the 
thermal atmospheric escape rates. We compare our results to 
the widely used energy-limited blow-off formula and discuss 
the relevance of our findings for "super-Earths" in general. 
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2 INPUT PARAMETER AND MODEL 
DESCRIPTION 

2.1 XUV fluxes of Kepler-11, GJ 1214 and 55 Cnc 

Unfortunately, the XUV fluxes of the host stars of the stud- 
ied "super-Earths" are not very weU constrained. The XUV 
emission of Kepler-11, a slightly evolved late-G star (Lis- 
sauer et al. 2011), is observationally unconstrained, because 
its large distance (> 600 pc) and high age (6-10 Gyr) pre- 
vent detection by current instruments. Hence, we estimate 
its XUV flux via its age by using power laws from Ribas 
et al. (2005), although we caution that application of these 
calibrations for stellar ages higher than about 6.7 Gyr (the 
highest age in their sample) might introduce additional un- 
certainties because of the extrapolation to an even older, 
slightly evolved G star. Adopting a mean age of 8 Gyr, this 
approach yields a logLxuv ~ 27.81 erg s~^, which corre- 
sponds to fluxes of Fxuv ~ 278, 204, 91, 61, 37ergcm-^ s"^ 
at the respective orbits of the Kepler-11 planets b-f. We es- 
timate the uncertainties of these values to be at least an 
order of magnitude due to the extrapolation of the calibra- 
tions and the just loosely constrained stellar age. 

Although GJ 1214 is located at a closer distance 
('^IS pc), the star has not been detected in X-rays up to now. 
This is probably due to the intrinsic faintness of this M4.5V 
star and its age of 3-10 Gyr (Charbonneau et al. 2009). 
Therefore, we use the relation between X-ray luminosity and 
age of Engle & Guinan (2011) calibrated for M dwarfs. This 
relation yields an X-ray luminosity of logLx ~ 26.8 erg s^^ 
for an adopted age of 6 Gyr. To correct for the unknown 
EUV emission, we assume that the contributions to the total 
XUV flux of both X-rays and EUV are ~50% each, typical 
for moderately active stars (e.g. Sanz-Forcada et al. 2011), 
but EUV can make up to 90% of the emission for inactive 
stars. This leads to an estimated XUV flux at GJ 1214b's 
orbit of Fxuv ~ 2200 erg cm~^ s~^. As for Kepler-11, the 
uncertainty of this value also corresponds at least an order 
of magnitude because of the large range of possible stellar 
ages, the intrinsic spread of about one order of magnitude of 
coeval stars of the same mass, and the very crude correction 
applied for the unknown EUV emission. 

For 55 Cnc, a nearby (~12.5pc) G91V star (von Braun 
et al. 2011), X-ray emission was detected by ESA's XMM- 
Newton space observatory (Sanz-Forcada et al. 2011). These 
authors used X-ray spectra of numerous planetary host 
stars to extrapolate their EUV emission by using calibra- 
tions of stars with well-determined emission measure dis- 
tributions (EMD). They estimated a total XUV luminos- 
ity within the wavelength region 5-920A of log 
27.551q'42 ergs~^, which translates to an XUV flux of 
Fxw ~ 4913 erg cm~^ s~^ at the orbit of 55 Cnc e. The 
uncertainty is 50%, mainly due to the uncertainties in 
the extrapolated transition region EMD. Although X-rays 
might contribute somewhat to the heating of the lower ther- 
mosphere, from our analysis of the available data we assume 
that for the particular atmospheric escape of the investi- 
gated planets is mainly driven by the stellar EUV emission. 
This assumption is supported by the fact that for rather old 
stars like the ones studied here, however, the major portion 
of the XUV flux (80-90%) is emitted in the EUV range. 

The properties of the studied planets are summarized 
in Table 1, which also gives the XUV enhancement factors 



Ixvv, corresponding to the ratio of the adopted XUV fluxes 
at the planet's orbits normalized to the present XUV flux of 
the Sun at 1 AU (4.64 erg cm"^ s"^ Ribas et al. 2005). 



2.2 Energy absorption and model description 

The energy budget of upper planetary atmospheres is mainly 
governed by the heating of the bulk atmosphere due to the 
absorption of the XUV radiation (e.g. Bauer & Lammer 
2004) and in very close orbital distances also harder X-rays 
(Cecchi-Pestellini et al. 2009; Owen & Jackson 2012) by at- 
mospheric constituents, by heat transport due to conduction 
and convection and by heat loss due to emissions in the in- 
frared (IR). Because the XUV fluxes are based only rough 
estimates except for 55 Cnc e, and the atmospheric com- 
position of the studied planets is not very well constraint, 
we do not consider a wavelength dependent derivation of 
the heating function. But we note that for future studies of 
exoplanets with known XUV spectra and atmospheric com- 
position it might be important to investigate the influence 
of wavelength dependent absorption cross sections and their 
impact on the upper atmosphere structure. The volume heat 
production rate gxuv due to the absorption of the stellar ra- 
diation for given wavelengths and constituents can then be 
written as 

gxuv ~ rjnas,Jxvve~'^ , (1) 

where Jxuv is the energy flux related to the corresponding 
wavelength range outside the atmosphere, Ua is the appro- 
priate hydrogen absorption cross-section which lies typically 
within the range of ~ 10~^*-6 x 10~^* cm^) and r is the opti- 
cal depth in the upper atmosphere and n the number density. 
»7 is the so-called heating efflciency or fraction of absorbed 
stellar radiation which is transformed into thermal energy 
and lies within in the range of ~15-60 % (Chasseflere 1996a; 
1996b; Yelle 2004; Lammer et al. 2009; Leitzinger et al. 2011; 
Koskinen et al. 2012). 

The value of rj is also related to the availability of IR- 
cooling molecules such as CO2, H^, etc. and will then be 
closer to the lower value because less energy can be trans- 
ferred into heat due to cooling by IR-emitting molecules. In 
a recent study by Koskinen et al. (2012), who applied a ther- 
mosphere model that calculates the heating rate based on 
the absorption of stellar XUV radiation and photoelectron 
heating efficiencies and photochemistry to the hydrogen-rich 
gas giant HD 209458b, it was found that rj is most likely in 
the order of ~40-60 %. For hydrogen-rich exoplanets which 
are exposed to XUV fluxes JSlOO times higher compared to 
that of the present Sun most molecules in the thermosphere 
are dissociated so that IR-cooling becomes less important 
yielding 77 ~40-60 %. As this is the case for HD 208459b, 
IR-cooling by Hg" molecules can be neglected (Koskinen et 
al. 2007). However, for the Kepler-11 "super-Earths" which 
are exposed to XUV fluxes ~8-60 times higher than today's 
Sun, depending on the availability of potential IR-cooling 
molecules, rj may be closer to the lower value of 15%. 

As mentioned in Sect. 2, it is also necessary to check if 
escape by X-ray heating could be relevant for the studied 
"super-Earths". However, as X-ray heating usually dom- 
inates at young stellar ages (Owen & Jackson 2012), it 
is unlikely to be of relevance here because of the rather 
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high ages of the host stars. If we compare the X-ray lu- 
minosities of the host stars with the values necessary for 
having dominating X-ray escape (Owen & Jackson 2012; 
their Fig. 11), one can estimate that X-ray luminosities 
of at least 10^* or lO'^^'ergs"^ would be necessary for the 
close-in GJ 1214b and 55 Cnc e "super-Earths" and the 
further separated Kepler- llb-f planets, respectively, if we 
adopt the values for a Neptune-mass planet with a density 
of lgcm~^ as a lower limit. The X-ray luminosities of all 
studied "super-Earth" host stars are all well below these 
values (< lO^^ergs"^). Therefore, we assume that X-ray 
induced evaporation can be neglected or may play a mi- 
nor role for the studied "super-Earths" , and XUV radiation 
should be responsible for the main heating processes in their 
hydrogen-rich thermospheres. 

By averaging the XUV volume heating rate over the 
planet's dayside one obtains 



Jxuv(r,e)sm(e)de, (2) 



with the polar angle O and Jxvv{r, O) — Jxuve"^. Because 
the thermal escape rate is determined by the total energy 
deposited in the upper atmosphere, which is given by the 
integral of the XUV volume heating rate over the radial dis- 
tance within the simulation domain, the shape of gxuv has 
a minor effect to the atmospheric loss rate. Besides the un- 
certainties in the availability of possible thermospheric IR- 
coolers, other parametrical uncertainties such as the chosen 
XUV flux, a possible X-ray heating contribution, the stellar 
age, etc. are most likely within the investigated limits of the 
assumed 15-40 % r] range. 

For studying the response of the hydrogen-dominated 
upper atmospheres of our "super-Earths" to the stellar XUV 
flux we apply a non-stationary ID hydrodynamic upper at- 
mosphere model which is a further developed model based 
on Penz et al. (2008), and solves the system of the 1-D fluid 
equations for mass, momentum, and energy conservation in 
spherical coordinates. 

dpr^ 



dt 



+ 



dpvr 
dr 



0, 



dpvr^ , d[r\pv' + P)] 

H = pgr + zPr, 



dt 

p^f. J P_ 



dr 



(7-1) 



dvr 



dt 



with pressure 
P = -^kT, 



+ 



(7-1) 



dr 

pvr^g + qiivvr^, 



and gravitational acceleration, 
5 = "V$, 



G 



Mel 



(d- 



G 



2d3 



Mstar + Mpl 



(3) 
(4) 

(5) 
(6) 
(7) 

(8) 



Here, r corresponds to the radial distance from the plan- 
etary center, p, v, P and T are the mass density, radial 



velocity, pressure and temperature of the bulk atmosphere, 
respectively, mn is the mass of the hydrogen atoms, k is the 
Boltzmann constant, G is Newton's gravitational constant, 
and 7 is the polytropic index or the ratio of the specific 
heats. 

For exoplanets at very close orbital distances one cannot 
neglect gravitational effects which are related to the Roche 
lobe (Erkaev et al. 2007; Penz et al. 2008; Lammer et al. 
2009a). Although, most of our studied "super-Earths" have 
orbital locations closer than 0.1 AU, due to their low gravity 
compared to Jovian-type planets their upper atmospheres 
may expand to several planetary radii and probably reach 
the Roche lobe. Therefore, we include in our study also the 
Roche lobe related gravitational force, which we refer to as 
Roche lobe effects with planetary mass Mpi is the and mass 
of the host star Mstar and orbital distance of the planet 
d. The upper atmospheres of planets can mainly be in two 
regimes. In the first, the thermosphere is in hydrostatic equi- 
librium where the bulk of the upper atmosphere gas below 
the exobase can be considered as hydrostatic. In the sec- 
ond, the upper atmosphere may hydrodynamically expand 
and the bulk atmospheric particles can escape efficiently as 
a result of high stellar XUV fluxes, and/or a weak plane- 
tary gravitational fleld (e.g., Lammer 2012). Under certain 
conditions and for certain planetary parameters the upper 
atmosphere can hydrodynamically expand, but the bulk ve- 
locity may not reach the escape velocity at the exobase level. 
In such cases one can expect a hydrodynamically expanding 
upper atmosphere but the loss rate results in a less efficient 
Jeans-type escape, but not in hydrodynamic blow-off (Tian 
et al. 2005b; Tian et al. 2008a; Tian et al. 2008b). If the 
mean thermal energy of the upper atmosphere gases at the 
exobase level exceeds their gravitational energy, blow-off oc- 
curs (e.g., Opik 1963). 

Our simulation domain is limited by a lower boundary 
with a temperature To, density no and a thermal velocity 



^^0 



( kTo 



(9) 



at the homopause distance or base of the thermosphere. 
To estimate the temperature To at the lower boundary one 
does not need to apply a full radiative transfer calculation 
because the temperature in this altitude region differs not 
greatly from the equilibrium or effective skin temperature 



(10) 



of the planet (e.g., Kasting & Pollack 1983; Bauer & Lammer 
2004; Tian et al. 2005b) . The factor ^ is 4 if a planet rotates 
fast and 2 for slowly rotating or tidally locked planets. Be- 
cause the studied "super-Earths" are within Mercury's orbit, 
it is likely that they rotate slower compared to the Earth so 
that ^ may be <4. ctb is the Stefan-Boltzmann constant, A 
the albedo and 5* the stellar fiux density which is the amount 
of the incoming electromagnetic radiation of the host star 
per unit area at the planet's orbit location. For the tempera- 
ture To at the lower boundary of the studied "super-Earths" 
we use the estimated effective skin temperatures Tcq given 
in Table 1. 

For planets close to or inside the habitable zone of a 
Sun-type G-star such as Venus or Earth this temperature is 
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Table 1. Planetary and stellar model input parameters for Kepler-llb-f (Lissauer et al. 2011; 
|http: / /kepler.nasa.gov/Mission/discoveries/| , GJ 1214b (Cliarbonneau et al. 2009) and 55 Cnc e (Demory et al. 2012; Endl et 
al. 2012; Gillon et al. 2012). /xuv corresponds to the stellar XUV flux at the planetary orbits normalized to the solar value at 1 AU. 



Exoplanet 


d [AU] 


-fxuv 




Mpi [Me] 


Ppl [g cm 3] 


To [K] 


Kepler-llb 


0.091 


60 


1.97 


4.3 


3.1 


~900 


Kepler- lie 


0.1 


45 


3.15 


13.5 


2.3 


~833 


Kepler- lid 


0.159 


20 


3.43 


6.1 


0.9 


~692 


Kepler- lie 


0.195 


13 


4.52 


8.4 


0.5 


~617 


Kepler- 1 If 


0.25 


8 


2.61 


2.3 


0.7 


-544 


GJ 1214b 


0.014 


470 


2.67 


6.55 


1.87 


~475 


55 Cnc e 


0.0156 


1060 


2.17 


8.37 


4.56 


-2360 




1 10 
r/rp, 

Figure 1. Calculated XUV volume heating rates gxuv for 55 
Cnc e (dashed line), GJ 1214b (dotted line), and the 5 Kepler-11 
"super-Earths" (solid lines) as a function of distance in planetary 
radii for rj = 15% (a) and ri = 40% (b), respectively. 



^230-250 K which is reached at an altitude of about 120 
km. 

Because of their low mean densities, in the case of GJ 
1214b (Bean et al. 2010; Miller-Ricci & Fortney 2010; Net- 
telmann et al. 2010; 2011) and 55 Cnc e (Demory et al. 
2012), it may be possible that these planets have a deep 
H2O ocean which is possibly surrounded by a steam atmo- 
sphere. Such a scenario is not expected for the Kepler-11 low 
density "super-Earths" , these planets are most likely domi- 
nated by dense H/He envelopes (Lissauer et al. 2011; Ikoma 
& Hori 2012). 

If a planet has a steam atmosphere the base of the ther- 
mosphere could also act as a kind of cold trap where the sat- 



uration H2O mixing ratio reaches its minimum value. If a 
planet is dominated by a hydrogen-rich gas envelope or has 
different parameters compared to the Earth, the homopause 
altitude could also be somewhat higher, but the exobase lo- 
cation or escape rate at the top of the thermosphere should 
not greatly differ. 

In the case of H2/He dominated upper atmospheres 
there are three processes which can dissociate II2 molecules 
into hydrogen atoms. The first process is thermal dissoci- 
ation which becomes dominant when the atmospheric tem- 
perature reaches values >2000 K (Yelle 2004; Koskinen et al. 
2010). The second process is photo-dissociation of hydrogen 
molecules which becomes dominant when the stellar EUV 
flux reaches values which are >25 times the present solar 
value (Koskinen et al. 2010). The third possibility for the 
production of atomic hydrogen are photochemical reactions 
between CH4 and other hydrocarbons taking place below the 
homopause level (e.g., Atreya et al. 1986). A combination of 
photo-dissociation and photochemical reactions is responsi- 
ble that H atoms are the dominate species in the upper at- 
mosphere of Jupiter, Saturn, Uranus and Neptune. Because 
the studied "super-Earths" are much hotter compared to the 
hydrogen-rich Solar System gas and ice giants, are therefore 
exposed to higher photon fluxes and may also contain hy- 
drocarbons in their lower atmospheres (Kuchner 2003) we 
assume that atomic hydrogen is the dominant species above 
their homopause levels. 

In case the planet is surrounded by a dense hot steam 
atmosphere, the stellar XUV radiation dissociates the H2O 
molecules above the mesopause via 

H20 + h!^^0H + H (11) 

and hydrogen atoms should also be the dominant species in 
the upper atmosphere (Kasting & Pollack, Tian et al. 2005; 
Lammer 2012). For the number densities no we assume a 
hydrogen density of 10^^ cm~^ similar as calculated for Sat- 
urn's, Uranus' and Neptune's homopause levels (Yamanaka 
1995; Atreya 1999). This value is also comparable with the 
homopause hydrogen density of 5 x 10^^ cm~^ which was 
adopted by Tian et al. (2005b) in a study of a hydrogen- 
rich early Earth, but slightly lower compared to that of 
Jupiter. This density value corresponds to a H2O mixing 
ratio in Earth's atmosphere of 50 % which results, accord- 
ing to Kasting & Pollack (1983) and Tian et al. (2005b), also 
in an upper atmosphere which is dominated by hydrogen. 

The upper boundary of our simulation domain is cho- 
sen at 45iipi, but the results of our hydrodynamic model are 
considered as accurate only until the exobase level which is 
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Table 2. Calculated exobase parameters without Roche lobe effects and thermal escape and mass-loss rates for the five "super-Earths" 
of the Kcpler-fl system, GJ f2f4b and 55 Cnc e for heating efficiencies r) of 15% (upper part) and 40% (lower part). 



Parameter 


Kepler-llb 


Kepler-llc 


Kepler-lld 


Kepler-llc 


Keplcr-llf 


GJ 1214b 


55 Cnc e 


Blow-off 

^cxo / ^pl 
Texo [K] 

rtcxo [cm^^] 
Tth [s'M 
Mth [gs-i] 


yes 
~ 33 
~1170 
~ 6 X 10^ 
~ 5.6 X 10^1 
~ 9.5 X 10^ 


no 
~ 14 
~ 1000 
~ 8 X 10^ 
~ 1.5 X 10^1 
~ 2.43 X lO'' 


yes 
~ 27 
~ 415 
~ 5 X 10^ 
~ 6.0 X 10^1 
~ 1.0 X 10* 


yes 
~ 25 
~ 430 
4 X 10^ 
~ 6.5 X 10^1 
~ 1.1 X 10** 


yes 

~ 41 

~ 50 
~ 1.6 X lO"* 
~ 2.0 X 10^2 
~ 3.3 X 10* 


yes 
~ 35 
~ 4040 
~ 5 X 10^ 
~ 2.25 X 10^2 
~ 3.45 X 10* 


yes 
~ 136 
2230 
~ 1.75 X 10^ 
~ 9.5 X 10^2 
~ 1.6 X 10^ 


Blow-ofF 

"^exo / ^pl 
Texo [K] 

rtcxo [cm~3] 
Tth [s-l] 
A/th [g s-i] 


yes 
~ 40 
~1800 
~ 6 X 10^ 
~ 1.0 X 10^2 
~ 1.5 X 10** 


yes 
~ 16 

~ 1560 

~ 10* 
~ 4.5 X 10^1 
~ 7.5 X 10^ 


yes 
~ 40 
~ 1040 
~ 6 X 10^ 
~ 6.5 X 10-^1 
~ 1.0 X 10** 


yes 
~ 28 
~ 680 
- 5 X 10^ 
~ 1.2 X 10^2 
~ 2.0 X 10* 


yes 
~ 45 
- 55 
~ 2 X 10* 
~ 3.7 X 10^^ 
~ 6.35 X 10* 


yes 
~ 40 
~ 7125 
6 X 10^ 
~ 4.5 X 10^2 
~ 7.5 X 10* 


yes 
~ 168 
- 5530 
~ 1.6 X 10^ 
~ 1.7 X 10^3 
~ 2.8 X 10^ 



located where the mean free path Zc reaches the scale height 
H — (kTcico) / (mg) of the hydrodynamically expanding up- 
per atmosphere. Above the exobase level which separates 
the collision dominated thermosphere from the collision-less 
exosphere, hydrodynamics are not valid. If /? > 30 at the 
exobase level then an atmosphere is mainly bound to the 
planet and one can expect low thermal escape rates. Blow- 
off happens when l3 becomes <1.5 or the upper atmosphere 
reaches the Roche lobe distance (Opik 1963; Chamberlain 
1963; Bauer & Lammer 2004; Erkaev et al. 2007). For ^ 
values between ~1.5-30, Jeans escape occurs but the upper 
atmosphere can expand within the hydrodynamic regime. 

Because the classical Jeans formula is based on the 
isotropic Maxwellian distribution function, in case if the 
upper atmosphere expands dynamically but blow-off con- 
ditions are not established, we have the radial velocity of 
the outward flowing bulk atmosphere at the exobase, and 
thus a distribution function which is not isotropic. In such 
a case we calculate the Jeans-type escape rate by using a 
shifted Majcwellian function which is modified by the ra- 
dial velocity v obtained from the hydrodynamic code (e.g., 
Volkov et al. 2011). We show the dynamically expanding 
upper atmosphere structures only up to the exobase levels. 



3 RESULTS 

The modeled XUV volume heating rates gxuv for the seven 
"super-Earths" with heating efficiencies 77 = 15% and 40% 
are shown in Fig. 1. Depending on the planets orbit loca- 
tion, equilibrium or skin temperature, and gravity, the peak 
of the XUV flux is deposited close to the planet or at higher 
altitudes. Corresponding to the volume heating rates the 
expected atmosphere temperature, velocity and density pro- 
files as a function of planetary distance without (Table 2) 
and with the Roche lobe effects (Table 3) are then mod- 
eled. Table 2 shows the calculated exobase parameters rexo, 
Texo, H'exo, the thermal escape rates Lth and mass-loss rates 
Mth for the five "super-Earths" in the Kepler- 11 system, GJ 
1214b and 55 Cnc e for heating efficiencies of 15% and 40 %, 
respectively. For example, for Kepler-llf the exobase would 
be located at ~ 117i?® ~ 10i?jup ~ lf?0. In the case of 
GJ 1214b which orbits around a M4.5 dwarf star with a ra- 
dius of ~ O.21i?0 ~ 23i?0 the exobase would be located at 



~ 4.0iZGJi2i4. This comparison indicates that the inflated 
upper atmosphere of GJ 1214b would cover the whole star 
during a transit. In the case of 55 Cnc e the expanding bulk 
atmosphere would reach escape velocity and hence blow-off 
conditions far below the estimated exobase level, so that the 
exobase location in that particular case is not relevant for 
the estimation of the thermal escape rate. For a heating ef- 
ficiency ri of 15% blow-off would occur at ~ 20f?pi and for 
rj =40% at ~ 14iipi. Due to the lower gravity the combi- 
nation of XUV heating and thermal expansion the exobase 
levels of Kepler-llb, Kepler-lld, Kepler-lle, Kepler-llf, GJ 
1214b and 55 Cnc e move to distances which are far beyond 
the Li point tli- Therefore, one can not neglect Roche lobe 
effects for these particular "super-Earths" . The thermal at- 
mospheric escape rates for most planets are ~ 2 times higher 
if the heating efficiency is 40 %. For ri = 40% all "super- 
Earths" are in the blow-off regime. For the lower heating 
efficiency of 15%, only Kepler-llc is not in blow-off but 
experiences strong Jeans-type escape at the exobase level 
from its non-hydrostatic expanded upper atmosphere which 
is about a factor of 0.65 lower compared to the hydrody- 
namic outflow rate at the exobase level. The usual Jeans 
approach which assumes zero bulk velocity would yield a 
value which is a factor of 0.145 lower. 

Another interesting effect on the exobase temperature 
can be seen at the most distant "super-Earth", Kepler-llf 
where Texo cools adiabatically to a very low value of ~ 55 
K. However, as mentioned above it is important to note 
that the results in Table 2 neglect the Roche lobe and its 
escape-enhancing effects. If we consider this more realis- 
tic scenario and compare the exobase levels of the studied 
"super-Earths" with the location of the Li point tli shown 
in Table 3 one can see that all "super-Earths" except Kepler- 
llc experience escape due to Roche lobe overflow. This ge- 
ometrical blow-off is most extreme for 55 Cnc e and is also 
effective for GJ 1214b for which the distances of the respec- 
tive LI points rLi are about 4f?pi. Due to the low gravities, 
the enormous expansions of the upper atmospheres and the 
relatively close distances to their host stars, one cannot ne- 
glect the Roche lobe effects. Table 3 shows the calculated 
parameters, but with Roche lobe effects included at the dis- 
tance of the corresponding LI points. In such a case the 
atoms escape as they reach and overflow the Roche lobe. 

If we compare the escape rates which neglect the Roche 
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Table 3. Calculated upper atmosphere parameters at the Roche lobe location and thermal escape and mass-loss rates for the studied 
"super-Earths" for heating efficiencies tj of 15% (upper part) and 40% (lower part). Because the exobase level of Kepler-llc does not 
reach rLi, the planet's values at the exobase are given instead. 



Parameter 


Kepler-llb 


Kepler-llc 


Kepler- lid 


Kepler-llc 


Kepler-llf 


GJ 1214b 


55 Cnc e 


rhi/Rpi 


~ 17.95 


~ 18.05 


~ 20.24 


~ 20.96 


~ 30.24 


~ 4.26 


~ 3.82 


r-exo/i?pl 




~ 14 












Tli [K] 


-670 


Tcxo ~ 785 


~ 353 


- 325 


~ 270 


~ 1360 


~ 1415 


"■Ll [cm 3] 


~ 3.7 X 10* 


"cxo ~ 7.0 X 10^ 


~ 1.0 X 10* 


6.5 X 10^ 


~ 9.4 X 10* 


~ 1.8 X 10^ 


~ 1.0 X 10^ 


Ah [s-M 


~ 7.0 X 10^^ 


2.45 X 10^1 


~ 6.0 X 10^^ 


~ 6.5 X 10^1 


~ 2.5 X 10^2 


~ 4.5 X 10^2 


~ 9.8 X 10^2 


Mth [g s-'] 


~ 1.15 X 10* 


~ 4.0 X 10^ 


~ 1.0 X 10** 


~ 1.1 X 10* 


~ 4.0 X 10** 


7.5 X 10* 


~ 1.6 X 10^ 


Texo/Rpl 




~ 16 












Tli [K] 


~1190 


Texo ~ 1415 


~ 650 


~ 590 


~ 110 


~ 2520 


~ 2890 


"Ll [cm 3] 


~ 4.8 X 10* 


ricxo - 9.7 X 103 


~ 2.2 X 10* 


~ 1.2 X 10* 


~ 7.0 X 10* 


~ 1.5 X lO** 


8.8 X 10" 


ith [s-i] 


~ 1.26 X 10^2 


~ 7.35 X 10^1 


~ 1.5 X 10^2 


~ 1.5 X 10^2 


~ 2.7 X 10^2 


~ 9.7 X 10^2 


^ 1.4 X 10^3 


Mth [g s-l] 


~ 2.0 X 10** 


~ 1.3 X 10** 


~ 2.5 X 10** 


2.5 X 10** 


~ 4.5 X 10** 


1.5 X 10^ 


~ 2.3 X 10^ 



Table 4. Estimated hydrogen loss rates from Eq. (12) for energy-limited escape {rj = 100%) and heating efficiencies of 40% and 15%. 
Exoplanet rj = 100%: Lth [s"M V = 40%: Lth [s"^] V = 15%: ^th [s"^] 



Kepler-llb ~ 6.0 X 10^2 ^ 2.4 X 10^2 ^ g g x lO^^ 

Kepler-llc ~ 5.9 X 10^2 ~ 2.3 X 10^2 ^ 8.9 x lO^i 

Kepler-lld ~ 7.5 X 10=*2 ~ 3.0 x 10^2 ^ ^ ^ iq32 

Kepler-lle ~ 8.0 X 10^2 ^ 3 2 x 10^2 ^ 12 X 10^2 

Kepler-llf ~ 3.5 X 10^2 ^ ^ 4 x 10^2 ^ 5.3 x 10^^ 

GJ 1214b ~ 4.6 X 10^3 ~ I.9 x 10^^ ~ 7.0 x 10^2 

55 Cnc e ~ 1.7 X 10=** ~ 7.0 x lO^^ ^ 2.6 x 10^^ 



lobe effects shown in Table 2 with those of Table 3 where 
Roche lobe effects are included, one can see that the stel- 
lar tidal forces enhance the atmospheric escape. The escape 
rates are influenced more strongly at the close-in "super- 
Earths" GJ 1214b and 55 Cnc e compared to the Kepler-11 
planets where tli is further away from Vpi. One can also see 
that the temperatures are different at the exobase level for 
Kepler-llc without (Table 2) and with Roche lobe effects 
(Table 3). The reason is that with the Roche lobe effects, 
the radial flow velocity is larger, which results in a larger 
adiabatic cooling. One should also note that a number den- 
sity which is less or larger than the assumed lO^'^ cm~^ at 
the homopause level would result in a lower or higher escape 
rate. 

Figs. 2 and 3 have the Roche lobe effects included and 
show the calculated temperature, velocity and density pro- 
files as a function of planetary radii from the homopause 
distance up to the Roche lobe and in the case of Kepler- 
llc up to the exobase level. In Fig. 2 we consider a low 
heating efficiency rj of 15%, i.e. the availability of IR-cooling 
molecules in the upper atmosphere and in Fig. 3 we use a 
heating efficiency = 40% as expected for hydrogen-rich 
hot gas giants (e.g., Yelle 2004; Koskinen et al. 2012). 

One can see in Figs. 2a and 3a that the temperature 
structure of these "super-Earth" upper atmospheres is de- 
termined by a complex interplay between the XUV flux, the 
planet's mean density, and its skin or equilibrium temper- 
ature To ~ Teq. In the case of the hottest "super-Earth" 
55 Cnc e the maximum XUV flux is deposited at higher 
altitudes (see also Fig. 1), therefore, one can see that the 
temperature flrst decreases until the XUV radiation heats 
the atmosphere and expansion related to adiabatic cool- 



ing becomes relevant. For Kepler-llc which is the heaviest 
"super-Earth" with a cooler skin temperature the XUV flux 
is deposited at a closer distance so that the atmosphere is 
heated to a temperature peak which is above 2500 K. Then 
adiabatic cooling due to hydrodynamic expansion cools the 
upper atmosphere at large distances to temperatures below 
1500 K for an tj of 40% and <1000 K if 77 = 15%. One can 
also see that for the "super-Earth" which is furthest away 
from its host star, namely Kepler-llf, the temperature de- 
creases until the upper atmosphere reaches the Roche lobe 
distance. The XUV flux is too low to raise the temperature 
around the altitude where the maximum energy is deposited. 
Due to this different behavior, Kepler-llf is the only "super- 
Earth" in our sample where the energy-limited formula in 
Eq. (12) underestimates the escape rate. From Figs. 2b and 
3b one can see that the expanding hydrogen atoms have ve- 
locities of a few km s~^ but <10 km s~*. Figs. 2c and 3c 
show the corresponding number density proflles. 

After the investigation of the blow-off regimes and up- 
per atmosphere structures of the seven "super-Earths" we 
compare the calculated thermal escape rates with escape es- 
timates from the simple energy-limited escape formula (e.g. 
Watson, Donahue, & Walker 1981; Hunten, Pepin, & Walker 
1987), but introduce a heating efficiency i] and a potential 
energy reduction factor due to do the stellar tidal forces 
(Erkaev et al. 2007) 



K = 1 



+ 



2(rLi/i?pi) 2(rLi/i?pi)3' 
The equation can then be written as 
3?7/xuvJxuv 



Lth 



(12) 



(13) 
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Figure 2. Modeled temperature, velocity and number density 
profiles of 55 Cnc e, GJ 1214b, Kepler-llb-f from the lower ther- 
mosphere up to the Roche lobe distance rLi given in Table 3 for 
a low heating efficiency with r) = 15%. 



where /xuv is the XUV enhancement factor given in Table 
1, Fxvv is the present time solar XUV flux at 1 AU, G is 
Newton's gravitational constant, ppi is the planetary density 
and rriH is the mass of a hydrogen atom. The results are 
shown in Table 4 for the energy-limited case with rj = 100% 
and heating efficiencies of 40% and 15%. 

The Roche lobe potential energy reduction factor K 
plays only a role for GJ 1214b and 55 Cnc e because the 
Li point of the Roche lobe is much closer to the planetary 
surfaces as compared to the other "super-Earths". K is 1.53 
for GJ 1214b and 1.63 for 55 Cnc e and ^ 1 for all five 
Kepler-11 planets. If we compare the thermal escape rates 
estimated by the simple Eq. (12) with those obtained by the 
hydrodynamic model in Table 3, one can see that, besides 
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Figure 3. Modeled temperature, velocity and number density 
profiles of 55 Cnc e, GJ 1214b, Kepler-llb-f from the lower ther- 
mosphere up to the Roche lobe distance tli given in Table 3 for 
a heating efficiency T) = 40%. 



of Kepler- 1 If, the atmospheric escape rates are only slightly 
overestimated. A difference of a factor 2-3 is certainly within 
all the other uncertainties such as the XUV flux enhance- 
ment. For that reason, our study indicates that the simple 
formula given in Eq. (12) is useful for mass-loss estimates of 
hydrogen-rich exoplanets during evolutionary time scales, if 
the particular planet is exposed to a high XUV flux and/or 
is located at close orbital distance. However, Eq. (12) should 
not be used if the atmosphere of a planet is not within the 
blow-off regime. By comparing the exobase location with the 
distance of the Li point of the most massive "super-Earth" 
Kepler-llc, one can see that the exobase does not reach the 
Roche lobe distance. By considering the Roche lobe effects 
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Figure 4. Comparison of the critical temperature Tc without 
(dashed line) and Tc{x) with Roche lobe effects (dotted line) at 
Kepler-llc as a function of planetary radii. The symbol "+" corre- 
sponds to Texo related to the low and higher rj values. The vertical 
solid line corresponds to the Roche lobe distance for the Li point 
and the shaded area marks the exobase levels of Kepler-llc for 
heating efficiencies rj of 15% and 40%, respectively. 



one can see from Table 3 that the exobase temperature Tcxo 
of Kepler-llc is ~785 K for a heating efficiency rj of 15% 
and ~1415 K if is 40%. The classical critical temperature 
jTc for blow-off can be written as (e.g. Hunten 1973) 



2GmHM, 



pi 



(14) 



By using the corresponding parameters for Kepler-llc given 
in Tables 1 and 3, is ~1490 K for »7=15% and ~1350 K 
if 77=40%. The corresponding exobase temperatures Toxo for 
the particular planet are ~785 K and 1415 K, respectively. 
By using this criteria one finds that the upper atmosphere 
of Kepler-llc would be in the blow-off stage for the higher 
heating efficiency but not for the lower value of 15%. How- 
ever, Erkaev et al. (2007) showed that stellar tidal forces 
influence also the critical temperature Tc so that under cer- 
tain conditions blow-off can occur. If one follows the deriva- 
tions related to the potential energy difference between the 
exobase level and the Roche lobe, given in Erkaev et al. 
(2007) we obtain the following equation 



Tc{x) = Tjup 



with 



and 



(15) 



(16) 



Here Tjup ~ 1.45 x 10^ K is the critical temperature for 
the onset of blow-off at the radius of Jupiter (i.e. when 
^cxo = Rpi), Mjup and i?jup are the mass and radius of 
Jupiter. Fig. 4 compares the temperature profiles Tc and 
Tc{x) of Kepler-llc as a function of planetary radius with- 
out and with Roche lobe effects. From this one can see that 
Roche lobe effects reduce the critical temperature Tc{x) at 
the exobase level to ~80 K for a heating efficiency 77 of 15%, 
resulting in an environment where Kepler-llc's upper at- 
mosphere will experience blow-off. But is interesting to note 
that Eq. (12) overestimates in that particular case the at- 
mospheric escape rate by a factor of ~3.6. 

Compared to the mass-loss estimations of "hot 



Jupiters" such as HD 209458b which are in the order of 
~2-5xl0^'' g s"^ (e.g. Yelle 2004; Tian et al. 2005a; Garcia 
Muiioz 2007; Penz et al. 2008; Lammer et al. 2009a; Koski- 
nen et al. 2012) the mass-loss of the studied "super-Earths" 
is ~10-100 times lower. The main reasons are that the stud- 
ied "super-Earths" are exposed to lower XUV fiuxes and the 
mean density of HD 209458b is lower. 



4 CONCLUSION 

We studied the blow-off regimes of seven hydrogen-rich 
"super-Earths" with known sizes and masses by applying a 
ID hydrodynamic upper atmosphere model which includes 
tidal forces. The upper atmosphere temperature profiles of 
these hydrogen-rich "super-Earths" are determined by a 
complex interplay between the equilibrium or skin temper- 
ature of the particular planet, the homopause number den- 
sity, the stellar XUV flux, the height where the maximum 
XUV energy is deposited, adiabatic cooling the mean den- 
sity of the particular planet, and the Roche lobe distance. 
The upper atmospheres of Kepler-llb, Kepler-lld, Kepler- 
lle, Kepler-llf, GJ 1214b and 55 Cnc e expand up to the 
Li points of their Roche lobes, which enhances their atmo- 
spheric escape rates. Kepler-llc is the only "super-Earth" 
in our sample where its exobase level does not reach the 
Roche lobe. The thermal mass-loss rates of GJ 1214b and 
55 Cnc e are ~10 times lower and about 100 times lower 
for the Kepler- 11 "super-Earths", respectively, compared to 
a typical "hot Jupiter" such as HD 209458b. By compar- 
ing the escape rates obtained by the hydrodynamic upper 
atmosphere model with the modifled energy-limited equa- 
tion given in Erkaev et al. (2007) and Lammer et al. (2009) 
we found that besides of Kepler-llc and Kepler-llf and 55 
Cnc e the escape rates are overestimated by about a factor 
two, which lies within all other parametrical uncertainties. 
However, our study also indicate that the energy-limited 
equation should only be applied for hydrogen-rich exoplan- 
ets which have a hot skin or equilibrium temperature and/or 
are exposed to high XUV flux values. From our results one 
can expect that the mass of the investigated "super-Earths" 
will not be affected strongly by atmospheric mass-loss dur- 
ing their remaining lifetimes. 
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